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ABSTRACT 

One of the important microlensing applications to stellar atmospheres is the study of spots on stellar surface 
provided by the high resolution of caustic-crossing binary-lens events. In this paper, we investigate the charac- 
teristics of spot-induced perturbations in microlensing light curves and explain the physical background of the 
characteristics. We explore the variation of the spot-induced perturbations depending on various parameters 
characterizing the spot and investigate how well these parameters can be retrieved from observations in high- 
cadence future lensing surveys. From this, we find that although it would not be easy to precisely constrain the 
shape and the surface brightness contrast, the size and location of the spot on the stellar surface can be fairly 
well constrained from the analysis of lensing light curves. 
Subject headings: gravitational lensing 
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1. INTRODUCTION 

Over the past decade, microlensing has developed into a 
powerful tool to study stellar ast rophysics, esp ecially stellar 
atmospheres. See the review of IGouldl d200lb . Microlens- 
ing application to stellar atmosphere is possible for events 
where the source star crosses the caustic. The caustic rep- 
resents the source position at which the lensing magnifica- 
tion of a point source formally becomes infinity. Then, as 
the caustic passes over the face of the star, different parts 
of the star are strongly magnified at different times, making 
it possible to probe detailed structures on the surface of the 
source star. One of the stellar properties that can be investi- 
gated by the high resolution provided by caustic-c rossing mi- 
crolensing events is the limb darke ni ng of stars (|Witj [19951 



Bogdanov & Cherepashchukl [19951: iGaudi & Gouldl 1 19991: 
Albrow et all 119991 1200 U lAfonso et al.1 l2000t IGouldl BOO it 



Hevrovskv 2003; Fields et al. 2003; Abe et al. 2003; Dominik 
2004t[Kubas et aU l2005UCassan et al.ll2006l) . 

Another application of microlensing to stellar at- 
mosphere is the st udy of irregular surf ac e stru ctures 
such as spots. lHevrovskv& Sasselov (2000) and 
iHendrv. Brvce & Valls-Gabaudl (120021) investigated the 
sensitivity to spots for single-lens events and found that the 
spot signal can be detecte d for events whe re t he lens transits 
the fac e of the source star. lHan et all (120001) and Chang & Hani 
(2002) demonstrated the feasibility of spot detections from 
the observation of caustic-crossing binary-lens events. 

Despite the feasibility of spot detections demonstrated by 
theoretical studies, there exists no firm detection of a spot sig- 
nal for any of the microlensing events detected so far. This 
might be due to the small fraction of stars with spots, but it is 
more likely that the current lensing surveys are not sensitive 
enough to catch the signal. For a caustic-crossing binary-lens 
event, the duration of a spot-induced signal is 
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where l, is the fraction of the spot size relative to the source 
size, p* is the source radius normalized by the Einstein radius 

1 corresponding author 



9e, ?e is the Einstein time scale, and <fi is the angle between 
the source trajectory and the caustic line. Considering that 
?e ~ 20 days for a typical Galactic event, the time scale is 
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Then, the spot-induced perturbation lasts only a few hours 
even for an event associated with a giant source star. On the 
other hand, the observational cadence of the current lensing 
surveys is several times per night, which is far less than that 
required to detect the spot-induced perturbations. 

However, the situation will be greatly different in future 
lensing surveys that will continuously survey wide fields at 
a high cadence using very large format imaging cameras. The 
OGLE collaboration recently upgraded their camera to widen 
the field of view from 0.4 deg 2 to 1.4 deg 2 . The MOA col- 
laboration plans to upgrade the telescope to a field of view 
of 4 deg" (T. Sumi 2009, private communication). The 'Ko- 
rea Microlensing Telescope Network (KMTNet)' is a sur- 
vey exclusive for microlensing that plans to achieve — 10 
minute sample using a network of three 1.6 m telescopes 
each of which is equipped with a camera of a 4 deg 2 field 
of view. With these upgraded and new instruments, the next- 
generation surveys will be able to resolve the short time-scale 
perturbations induced by spots. 

In this paper, we further explore the feasibility of mi- 
crolensing studies of stellar spots. We investigate the char- 
acteristics of spot-induced perturbations in microlensing light 
curves and explain the physical background of the character- 
istics. We also investigate the feasibility of constraining the 
physical parameters of spots from the observations in future 
lensing surveys. 

2. SPOT-INDUCED PERTURBATIONS 

2.1. Characteristics 

Figure Q] shows the light curve of an example caustic- 
crossing binary-lens event occurring on a source star with a 
spot and the resulting pattern of the spot-induced perturba- 
tion. From the investigation of the perturbation, one finds 
the two following important characteristics of the perturbation 
pattern. First, there exist both positive and negative deviations 
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in lensing magnification from the light curve of the event oc- 
curring on an unspotted star. Second, the perturbation pattern 
is asymmetric with respect to the maximum deviation. 

For an event occurred on a spotted star, the lensing magni- 
fication is represented by 



A(x,y)dxdy — (1 — C s ) J s A(x,y)dxdy 



£*-£. 



(3) 



where A(x,y) represents the point-source magnification at a 
source position (x,y), C s is the flux ratio between the spot and 
unspotted regions, and and E. represent the surface areas 
of the star and the spot, respectively. Here we assume that the 
brightness distribution on the source star surface other than 
the spot region is uniform. 2 Then, the first term of the nu- 
merator of the right side of equation (f3]l represents the contri- 
bution of the unspotted region to the magnification (non-spot 
term), while the second term represents the contribution of the 
spot region (spot term). 

The negative deviation in the spot-induced perturbation oc- 
curs when the magnification of the spot region is very high. 
In this case, the decrease of the numerator in equation (01 by 
the spot term is more important than the decrease of the de- 
nominator by the spot area, resulting in a negative deviation. 
The magnification of the spot region is maximized when the 
spot is on the caustic. 

The positive deviation occurs when the caustic is over the 
source surface but the spot is away from the caustic. In this 
case, the magnification at the position of the spot is not high 
and the non-spot term in equation (f3]i dominates over the spot 
term. Then, the decrease in the numerator is negligible while 
the decrease of the denominator remains the same regardless 
of the spot position. With the increment of the denomina- 
tor, the magnification is higher than the magnification of an 
unspotted event, resulting in a positive deviation. 

The asymmetry of the perturbation pattern is caused by the 
difference in magnification patterns between the outside and 
the inside regions of the caustic. In the outer region, the 
magnification plummets as the distance d from the caustic in- 
creases, while the magnification inside the caustic decreases 
smoothly as A oc c/" 1 / 2 . Then, the spot term in equation (f3]i is 
important in a very narrow region outside of the caustic but 
in a much wider region inside the caustic. As a result, the 
deviation pattern is asymmetric. 

With the understanding of the physical background of the 
patterns of spot-induced perturbations, the locations of several 
important turning points in the perturbation can be explained. 
For the case of the example event presented in Figure Q] there 
exist three such points, each of which is marked by a number. 

1 . The position marked by "1" corresponds to the moment 
at which the source surface enters the caustic. Differ- 
ential magnification becomes important from this mo- 
ment. The highly magnified parts on the source star sur- 
face are unspotted regions and thus the non-spot terms 
dominates over the spot term, resulting in a positive de- 
viation. 

2. The position marked by "2" corresponds to the mo- 
ment of spot's entrance into the caustic. From this mo- 
ment, the spot term becomes important and the devia- 
tion rapidly drops into negative values. 

2 For the effect of limb darkening on the spot-induced perturbation, see 
Hendry, Bryce & Valls-GabauJ (20021) 
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FIG. 1 . — Light curve of a caustic-crossing binary-lens event occurred on 
a source star with a spot (middle panel) and the the residual from that of the 
event occurring on an unspotted source star (lower panel). The light curve is 
in lensing magnification and the residual is the difference in magnification. 
The upper panel shows the geometry of the lens system where the thick linear 
structure represents the caustic and the straight line with an arrow represents 
the source trajectory. The grey scale shows the point-source magnification 
as a function of position and brighter tone represents a higher magnification. 
The inset shows the wider view of the geometry. The circles with different 
colors represents the source at the times marked by numbers in the lower 
panel. The spot of each source is also marked. The three light curves in 
the middle panel represent those of events occurred on source stars with a 
spot (solid curve), without a spot (dashed curve), and a point source (dotted 
curve), respectively. The time is expressed with respect to the time of the 
caustic crossing of the source center, t cc , in units of the caustic-crossing time 
scale At cc . 

3. The negative deviation is maximized when the spot 
term is maximized and this corresponds to the position 
marked by "3". We note that the magnification inside 
and outside of the caustic is very different and thus the 
moment of maximum negative deviation occurs not at 
the moment when the center of the spot is exactly on the 
caustic but when the center is slightly shifted toward the 
inside of the caustic. 



2.2. Variation 

Based on the fundamental scheme of the deviation de- 
scribed in the previous subsection, spot-induced perturbations 
take different shapes depending on various factors. These fac- 
tors include the size, shape, location on the source surface, 
orientation, brightness contrast with respect to the unspotted 
region, etc. 

Figure [2] shows how the spot-induced perturbation depends 
on various factors. Under the simplified approximation of a 
single elliptical spot with a uniform surface brightness, we 
parameterize the spot size as the semi-major axis a s normal- 
ized by the source radius, the shape as the axis ratio R s of the 
ellipse, the location on the source as the polar coordinates (r s , 
# s ) of the center of the ellipse with respect to the source center, 
the orientation as the angle ip s between the semi-major axis of 
the spot and the source trajectory, and the brightness contrast 
as the flux ratio C s between the spot and unspotted regions. 

3. CHARACTERIZATION 
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FIG. 2. — Variation of the spot-induced perturbations. The top panel shows the geometry of the lens system with respect to the caustic produced by a binary 
lens system with a mass ratio q = 0.5 and a projected separation normalized by the Einstein radius between the components of s = 0.7. The source radius relative 
to the Einstein radius is p* = 0.02. The subsequent panels show the variation depending on the spot size, the shape, the location, the orientation, and the contrast 
in brightness between the spot and unspotted regions, respectively. 



In the previous section, we explored the variation of the 
pattern of spot-induced perturbations depending on various 
parameters that characterize the spot. In this section, we in- 
vestigate how well these parameters can be retrieved from the 
analysis of light curves of event to be observed in future lens- 
ing surveys. For this, we estimate the uncertainties of the spot 
parameters by fitting an example caustic-crossing binary-lens 
event produced by simulation under the observational condi- 
tions of planned future lensing surveys. 

We choose the lensing parameters of the example event by 
adopting those of a typical Galactic bulge event occurred on 
a giant star. The adopted values of the Einstein time scale is 
?e = 20 days. The projected separation normalized by the Ein- 
stein radius and the mass ratio between the lens components 
are s = 0.7 and q = 0.5, respectively. The apparent baseline 
magnitude of the source is / = 15 and the fraction of the flux 
from the lensed star is 70% and the rest comes from blended 
stars. The source trajectory with respect to the caustic is 
same as shown in the top panel of Figure [2] The values of 
the parameters characterizing the spot properties are a s = 0.3, 
R s = 0.5, ip s = 45°, r s = 0.5, 9 S = 45°, and C s = 0.5. 

The event is assumed to be continuously observed in the / 
passband by using a network of telescopes with a cadence of 



6 times per hour. The exposure of each observation is 30 sec- 
onds. The aperture of each telescope is 1.6 m and the quantum 
efficiency of the detector is 0.8. Given the specification of the 
instrument and the magnitude of the source and blend, we set 
the photometric uncertainty by assuming that the photome- 
try follows photon statistics with a 1 % systematic uncertainty. 
We also assume that the photometry is Gaussian distributed. 
Figure|3]shows the caustic-crossing part of the light curve pro- 
duced by the simulation (upper panel) and the residual from 
the unspotted light curve (lower panel). 

The uncertainties of the spot parameters are estimated by 
fitting the light curve produced by simulation. In the fitting 
process, we search for a set of parameters describing the per- 
turbation by minimizing x 2 m tne parameter space. We use 
a Markov Chain Monte Carlo method for the \ 2 minimiza- 
tion. Since a spot-induced perturbation occurs in a small re- 
gion during the caustic crossing, its effect on the other lens- 
ing parameters is negligible. We, therefore, set the lensing 
parameters fixed in the fitting process. 3 We hold two of the 

3 Another reason for the use of this approximation is the limitation of com- 
putation time. Modeling binary-lensing light curves requires to include many 
parameters. To describe spot-induced perturbations, at least 5 additional pa- 
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FIG. 3. — Light curve of the event produced by simulation that is used for 
the estimation of the uncertainties of the spot parameters. We note that the 
light curve is expressed in magnitude. The time ?o represents the moment of 
the closest approach of the source star to the center of mass of the binary lens. 



TABLE 1 

Uncertainties of Spot Parameters 



property 


parameter 




uncertainty 


semi-major axis 


a s 




= 0.04 


axis ratio 




&R S 


= 0.2 


orientation 


Vs 


°"</>s 


= 30° 


location of the source surface 






,<r e ,) = (0.02,5°) 


surface-brightness contrast 


c s 


CT Cs 


= 04 



spot parameters fixed at a grid of values, while the remaining 
parameters are allowed to vary so that the model light curve 
results in a minimum \ 2 at eacn grid point. We set a s and 
R s as the grid parameters because they characterize the two 
most important properties of a spot. Another important rea- 
son for using grid parameters is that it allows to use multiple 
CPUs by allocating computations of different sections of the 
grid-parameter space different CUPs. Once a series of mod- 
els for the individual sets of the grid parameters are obtained, 
we estimate the uncertainty of each parameter from the distri- 
bution of the parameter obtained by giving weight to the % 2 
difference from the value of the best-fit model. 

For the computation of lensing magnification affected 
by spot perturbation, we use the ray-shooting technique 
(ISchneider & Weissl fl986l iKayser. Refsdal. & StabeiH 119861: 
Wambsgan ssl 19971) . In this technique, a large number of light 
rays are uniformly shot from the observer plane through the 
lens plane and then collected on the source plane. We ac- 
celerate this process by restricting the region of ray-shooting 
only in the region around the caustic and using a simple semi- 
analytic approximation in other parts of the source plane. In 
addition, we keep the information of the positions of the light 
rays arriving at the target in the buffer memory of the com- 
puter so that it can be readily used for fast computation of the 
magnification. 

We present the distributions of A\ 2 for the individual pa- 
rameters are needed. As a result, it is difficult to fir light curves letting all 
parameters vary. Fortunately, the variation of the lensing parameters on the 
determined spot parameters will be very minor and thus will not affect the 
result of the analysis. 
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FIG. 4. — Distributions of Ax 2 for the spot parameters. The regions marked 
in red, yellow, and green regions represent those of A\ 2 < 1, 4, and 9, re- 
spectively. 

rameters in Figure [4] and list the determined uncertainties in 
TableQ] From the uncertainties, we find the following factors. 
First, the size of the spot can be well constrained with uncer- 
tainties of cr fls ~ 0.05. We judge the reason for this is due to 
the existence of the discrete structures on the spot perturba- 
tion such as the turning points mentioned in section 2. These 
structures depends on the spot size and thus help to constrain 
the size. Second, the location of the spot on the surface of 
the source star is also well constrained. This is because the 
source surface is swept by the caustic two times, i.e. at the 
entrance and exit. For the individual caustic crossings, the 
incidence angles of the source respect to the caustic are dif- 
ferent and thus the location of the spot is well constrained. 
On the other hand, the shape of the spot is expected to be de- 
termined with a relatively large uncertainty. Same is true for 
the surface-brightness contrast. This can be understood from 
the similarity in the variation of the perturbation pattern de- 
pending on the two spot parameters of R s and C s as shown in 
Figure |2] where it is found that the variation of the perturba- 
tion with the increase of R s is imitated by the variation with 
the increase of C s . This implies that an observed perturba- 
tion suffers from the degeneracy between the two parameters 
and thus results in large uncertainties of both parameters. The 
expected uncertainty of the orientation angle of the spot is 
~ 10°. 
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The uncertainties estimated in this work might be underesti- 
mated compared to the actual values in the future experiments. 
Possible causes of the underestimation would be the simpli- 
fication of the spot as a uniform ellipse and the disregard of 
detailed structures such as umbra and pen-umbra. Consid- 
ering that spots would occupy a small fraction of the source 
surface and the estimated uncertainty of the spot-size is ~ 5% 
of the stellar size, we think that the additional uncertainties 
caused by the disregard of very detailed structures would not 
be important. Another cause would be the disregard of the 
stellar rotation. Due to the stellar rotation, the locations of 
the spot on the stellar surface at the times of the caustic en- 
trance and exit are different, causing additional uncertainties, 
especially in r s and 6 S . We note, however, that most of the 
events for which spots can be detected will be events asso- 
ciated with late-type giant source stars for which the angular 
rotation speed is small. Therefore, the uncertainty would not 
be seriously different from our estimation. 

4. CONCLUSION 

We investigated the characteristics of spot-induced pertur- 
bations in microlensing light curves and explained the phys- 



ical background of the characteristics. We also explored the 
variation of the spot-induced perturbations depending on var- 
ious parameters characterizing the spot and investigated how 
well these parameters can be retrieved from observations in 
high-cadence future lensing surveys. From this, we found 
that the size and location of the spot on the stellar surface can 
be fairly well constrained from the analysis of lensing light 
curves, although it would not be easy to precisely constrain 
the shape and the surface brightness contrast 
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